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Information  can  be  manipulated  and  processed  by  a  nonlinear  dynamical  system.  This 
work  investigates  the  properties  of  such  nonlinear  optical  system's  based  up'on  real-time 
holographic  materials  know'n  as  photorefracti  ve  media.  Our  goal  is  to  understand  how 
and  when  a  dynamical  system  can  self-organize,  or  learn  on  its  own  to  process 
information  in  a  useful  way.  We  have  demonstrated  systems,  for  example,  that  can 
separate  multiple  fiber  optic  communication  channels  or  learn  to  classify  incomini: 
images.  We  have  discovered  that  such  self-organizing  systems  will  exhibit  a  phase" 
transition  with  the  information  content  of  the  input  as  the  critical  parameter.  The  phase- 
transitions  are  analogous  to  those  common  to  thermodynamic  systems.  In  order  to  better 
provide  an  analytical  basis  of  treatment  of  our  systems  we  have  developed  various 
analytical  tools,  in  particular  .stability  analysis  tools.  We  have  also  looked  at  the 
photorefractive  version  of  well  known  nonlinear  dynamical  problems  such  as  the 
formation  of  hexagonal  patterns  in  pumped  photorefractive  media. 
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I.  Overview 


Research  for  the  period  June  1992  to  June  t9Q3  ha-:  encompassed  sevciai  topics  in  coir.petiti  .c 
dynamics  and  self  organization  in  photorefractive  systems; 

1)  Stability  analysis  of  photorefractive  circuits. 

2)  Spatio-temporal  dynamics  in  a  photorefractive  ring  oscillator. 

3)  Transverse  instability  of  counterpropagating  beams  in  photorefractive  media. 

A  short  synopsis  of  results  in  each  of  these  areas  is  included  in  the  rest  of  this  section.  v.hile 
additional  details  are  given  separately  for  each  topic  in  the  following  sections. 

It  was  demonstrated  over  the  past  several  years  that  photorefractive  ring  resonators  and  circuits  of 
coupled  photorefractive  resonators  are  suitable  for  implementing  a  wide  variety  of  competitiv  e  and 
cooperative  dynamics  [1-7]  including  specifically,  the  bistable  ring  resonator  w  ith  gain  and  loss 
[1,2],  "winner-takes-all”  and  "voting  paradox  dynamics"  [3],  a  photorefractive  flip-flop  [4],  a 
ring-resonator  with  history  dependent  feedback  [5],  and  self-organizing  feature  extractors  that 
operate  on  frequency  or  time  multiplexed  spatial  patterns  [6,7].  These  systems  were  demonstrated 
experimentally  before  a  full  theoretical  analysis  was  completed.  We  have  now  completed  stability 
analyses  of  the  flip-flop  [4],  and  the  self-organizing  feature  extractor  [6.7]  that  confirm  the 
experimental  observations.  A  detailed  analysis  of  the  flip-flop  [27]  clarifies  the  role  of  the  time 
constants  in  the  gain  and  loss  media  on  the  .system  stability.  Analysis  of  a  simplified  model  of  the 
feature  extractor  [8]  shows  that  the  desired  stationary  state  corresponding  to  extraction  of  the  input 
patterns  is  stable,  and  all  other  stationary  .states  are  unstable.  In  the  case  of  the  feature  extractor 
numerical  analysis  has  been  used  to  study  the  time  dependence  of  the  transition  to  a  .self-organized 
state,  and  the  role  of  noise  induced  symmetry  breaking  in  this  system. 

In  order  to  study  more  complex  self-organizing  systems  than  the  feature  extractor  we  hav  c  moved 
to  an  investigation  of  a  single,  transversely  continuous  ring  resonator  [9].  This  resonator  allows 
for  the  self-organized  formation  of  transversely  continuous  spatial  patterns,  whereas  in  the  coupled 
ring  resonators  the  complexity  of  the  self-organized  pattern  is  limited  by  the  number  of  coupled 
rings.  Wc  choose  to  work  with  an  imaging  resonator  that  admits  an  arbitrary  transverse  mode  as  an 
eigenmode.  In  this  situation,  where  the  passive  optics  play  only  a  minor  role  in  determining  the 


mode  structure  the  impact  of  the  nonlinear  elements  comes  to  the  forefront.  Our  initial 
experimental  results  demonstrated  how  photorefractive  gain  and  loss,  placed  in  Fourier  conjugate 
resonator  planes,  induce  the  continuum  of  transverse  modes  to  collapse  into  a  single  localized 
mode  of  near  Gaussian  form.  Since  the  imaging  ring  has  no  optical  axis  the  location  of  the 
localized  mode  is  arbitrar>',  and  is  seen  to  form  at  any  point  in  the  transverse  plane  of  the  resonator. 
In  fact,  due  to  imperfect  alignment  of  the  resonator,  the  localized  mode,  or  spot,  is  not  stationary, 
but  rather  drifts  in  the  transverse  plane  (10).  This  phenomena  appears  to  be  unique  in  the  field  of 
nonlinear  optics.  Analytical  and  numerical  studies  of  the  resonator  show  that  the  rate  of  spot 
motion  is  directly  dependent  on  the  amount  of  cavity  misalignment. 

Part  of  the  motivation  for  studying  nonlinear  dynamics  in  photorefractive  systems  is  the  possibility 
of  implementing  neuromorphic  information  processing  systems.  In  this  context  the  localized 
excitation  in  the  photorefractive  ring  resonator  represents  the  response  of  the  information 
processing  system  to  an  input  pattern  (the  resonator  pump  beam).  We  have  investigated  the 
response  of  an  idealized  resonator  w  ith  no  cavity  misalignment  to  time  multiplexed  input  patterns. 
We  find  from  numerical  simulations  that  input  patterns  with  a  large  degree  of  similarity  lead  to 
resonator  excitations  lying  close  together,  w'hile  dissimilar  input  patterns  lead  to  resonator 
excitations  spaced  far  apart.  It  appears  therefore  possible  to  use  this  resonator  concept  to  implement 
a  class  of  neural  network  algorithms  due  to  Kohonen.  known  as  topologv  preserving  mappings 
[111- 

While  the  motion  of  the  localized  excitation  is  interesting  in  it.>  own  right  as  an  example  of  novel 
dynamical  behavior,  it  is  undesirable  for  implementing  e.g.  a  topology  pre.serving  mapping  that 
should  result  in  a  stationary  output  when  presented  with  a  .statistically  stationary  input.  We  have 
therefore  investigated  various  ways  in  which  additional  nonlinearities  may  be  used  to  inhibit  the 
motion,  while  still  allowing  the  mode  to  form  at  an  arbitrary  position.  Most  of  this  work  has  been 
performed  numerically  using  a  3-D  simulation  program  that  allows  a  wide  variety  of  resonator 
configurations  to  be  investigated.  In  early  June  1993  we  succeeded  in  finding  a  configuration, 
requiring  one  additional  photorefractive  crystal,  that  stabilizes  the  spot  motion  on  the  computer, 
despite  cavity  misalignments.  Work  is  now  under  way  to  demonstrate  this  experimentally. 

In  addition  to  the  work  on  photorefractive  ring  resonators  we  have  studied  a  more  traditional  topic 
in  nonlinear  optics:  the  transverse  modulational  instability  of  counterpropagating  waves.  The 
modulational  instability  is  well  known  in  nonlinear  optics,  dating  back  to  some  Soviet  work  from 
the  late  1970's  [12].  The  nonlinear  stage  of  {he  instability  re.sults  in  the  formation  of  transverse 
patterns,  including  rings,  squares  and  hexagons,  and  has  been  observed  experimentally  in  atomic 
vapors  [13-15],  liquid  crystals  [16,17].  and  quite  recently  in  photorefractive  media  [18].  The 
intrinsically  slow  dynamics  of  photorefractive  media  wnich  simplifies  time  resolved  measurements, 
and  their  low  thresholds,  which  make  the  strongly  nonlinear  stage  of  the  instability  accessible, 
render  them  well  suited  to  further  experimental  investigation  of  these  phenomena.  The  likelihood 
that  photorefractive  media  will  come  to  contribute  to  better  understanding  of  the  nonlinear  aspects 
of  pattern  formation  in  this  system  motivates  the  theoretical  analysis  of  the  instability  threshold  that 
we  have  undertaken  as  a  first  step  towards  analyzing  the  nonlinear  stage  of  the  instability. 


The  analysis  is  complicated  in  photorefractive  media  by  the  fact  that  the  coupling  constant  is  in  n/ 
gentidl  cuiiijdcA.  whereas  it  is  purely  real  in  the  instantaneous  Kerr  media  that  have  been  the  -  y  Codes 
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subject  of  previous  theoretical  work.  A  paper  describing  the  instability  threshold  condition  due  to 
the  formation  of  transmission  gratings  in  photorefractive  media  has  been  accepted  for  publication 
[19].  A  second  theoretical  contribution  in  this  area  that  considers  the  reflection  grating  case  which 
was  experimentally  realized  [18],  is  nearing  completion  [20]. 


II.  Stability  analysis  of  photorefractive  circuits 


The  analytical  tools  necessary  for  detailed  analysis  of  the  stability  of  photorefractive  circuits  based 
on  two-beam  coupling  can  be  found  as  an  extension  of  the  methods  that  were  previously  developed 
for  treating  a  number  of  four-wave  mixing  configurations  in  photorefractive  media  [28].  The 
conceptual  basis  of  the  method  is  quite  straightforward:  perform  a  linear  stability  analysis  about  all 
possible  stationary  solutions  of  the  circuit  under  consideration,  in  order  to  isolate  the  stable 
solutions.  The  mechanics  of  carrying  this  program  out  becomes  involved  as  the  complexity  of  the 
circuits  increases.  The  analysis  is  lengthy  because  the  steady  state  solutions  almost  always  involve 
nonlinear  variations  of  the  optical  fields  and  photorefractive  index  gratings.  Thus  the  linear 
perturbation  analysis  must  be  performed  about  a  nonlinear  steady  state. 

The  simplest  configuration  we  wish  to  analyze  is  the  dynamics  of  two  beam  coupling,  which  is 
described  by  the  equations 


(11.1) 


where  s  and  p  are  the  amplitudes  of  the  signal  and  pumping  beams  respectively,  j?  is  the  amplitude 
of  the  grating,  T  is  the  coupling  constant,  T  «  is  the  characteristic  time  constant  of  the  medium 
and  Ij  =  |.v|^  +  IpI"  is  the  sum  of  intensities  of  the  interacting  beams.  The  well  known  steady  state 
solutions  of  Eqs.(II.l)  for  real  coupling  may  be  written  in  the  form 

KIL  =  IpoL  =kl,>exp(-r), 

^  1  +  ,  (II.2) 

M  = - , 

r +  exp(-r) 

where  r  =  j.v,) is  the  input  beam  intensity  ratio.  Equations  (11.2)  allow  for  any  sign  of  f  . 
Below  we  will  use  the  notation  T  s  r^;,  M  =  G  for  T  >  0  (when  the  signal  beam  experiences 
gain)  and  T  =  -F,.  M  ^  L  kn  F  <  0  (.when  ii  expeiiences  loss^ 
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Arbitrary  complex  perturbations  about  the  steady  state  solution  can  be  separated  into  purely  real 
and  purely  imaginary  contributions  and  their  evolution  can  be  described  independently.  Below  \kc 
restrict  ourselves  to  consideration  of  purely  real  perturbations,  since  they  turn  out  to  be  more 
important  foi  the  analysis  of  the  circuits  considered  here.  The  time  dependent  amplitudes  are 
therefore  written  in  the  form 


.s'(.v./)  =  .Sf,(.v)  +  Re(5.v(.v)exp(/if)). 
p{  .V.  ? )  =  p,)  ( .V )  +  Re  ( Sp(  .V )  e  xp(  fi )). 


(II. 


where  f  is  a  complex  frequency. 

Linearizing  Eqs.(II.l)  with  respect  to  5s  and  5p  around  the  stationary  solution  (11.2)  and  sohing 
the  resulting  equations  results  in  a  transmission  matrix,  that  describes  propagation  of  perturbations 
through  a  medium  by  coupling  their  input  and  output  amplitudes 


tll.4) 


where  the  2  x  2  transmission  matrix  elements  are  given  by  the  relations; 


T-(-)  exp(r/2) 

■'ll  ~ 
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A  =  1  +  rexp(r). 
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Figure  II.  1  The  photorefractive  flip-flop 


There  is  also  a  matrix  T  "  describing  the  propagation  of  perturbations  when  a  pump  beam  interact'' 
with  two  spatially  uncorrelated  signal  fields. 

With  the  help  of  these  matrices  the  stability  of  very  general  photorefractive  circuits  may  be 
analyzed.  The  details  turn  out  to  be  lengthy  and  only  a  summary  of  the  physical  implications  is 
given  below. 

The  photorefractive  flip-flop  is  shown  in  Fig.  II.  1  [4].  Conditions  for  existence  of  spatial 
bistability  based  on  steady  state  considerations  were  given  in  [4].  This  analysis  also  shows  the 
effect  of  varying  the  crystal  time  constants.  The  flip-flop  has  six  steady  states:  both  rings  off.  one 
ring  on  and  the  other  ring  off,  both  rings  on  equally  and  both  rings  on  with  unequal  intensities.  It 
turns  out  that  all  the  stationary  states,  including  the  flip-flop  state,  are  unstable  when  T,  »  r,.,  . 
So  it  is  necessary  to  have  a  fast  loss  crystal,  for  the  flip-flop  to  be  stable.  The  same  type  of 
constraint  applies  in  e.g.,  the  bistable  ring  resonator  [1,2]  where  the  off  state  is  unstable  w  hen  the 
loss  mechanism  is  too  slow. 


A  slightly  different  configuration  where  the  two  rings  share  a  common  gain  volume  has  also  been 
analyzed.  The  shared  volume  leads  to  additional  competition  between  the  rings  and  the  flip-flop 
state  is  now  stable,  while  all  other  states  are  unstable,  when  r,  »  r^,-.  The  experimental  realization 
of  the  flip-flop  [4]  resulted  in  partial  sharing  of  the  gain  volume  between  the  rings,  and  had 


A  simplified  model  of  the  self-organizing  feature  extractor  based  on  two  coupled  single-mode 
photorefractive  ring  oscillators  is  showm  in  Fig.  11.2.  The  device  that  was  experimentally 
demonstrated  was  based  on  two  coupled  multi-mode  resonators,  each  containing  an  additional 
photorefractive  crystal  that  ensured  that  the  multi-mode  ring  behaved  as  though  it  were  a  single 
mode  ring.  In  order  to  make  the  problem  analytically  tractable  each  multi-mode  ring  w'ith  its  extra 
photorefractive  crystal  i:,  ;r;ndcled  as  an  empty  single  mode  ring.  The  .set  of  equations 


fX,  I  , 


()fy 


=-i.’ 


1’  .V 


'  ik  ' 


T  =  TJI,.  /7=SK+/'<0 

*; 

captures  the  essential  dynamics  of  this  system.  In  (11.6)  5,*  is  the  field  amplitude  of 


Output  Output 
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Figure  II. 2  Simplified  model  of  a  photorefractive  feature  extractor  consisting  of  two  single-mode 
resonators  with  shared  gain  volume. 

signal  k  in  ring  /.  is  the  field  amplitude  of  the  input  pump  beam  in  spatial  mode  j  and  temporal 
mode  k,  is  the  refractive  index  grating  coupling  ring  i  to  pump  spatial  mode  k.  r  is  the 
photorefractive  time  constant  and  is  the  total  optical  intensity. 

These  equations  are  supplemented  by  the  boundary  conditions 

/?,,(- =  0,0  = P22(z  =  0,t)  =  Po22^  p,2(2  =  0,r)  =  /72i(t:  =  0,0  =  0 

5„(O,/)  =  e-^'/^5„(/,0,  .r„(0,0  =  e-^''^Vu(/T),  g,U,0)  =  0.  (II.7) 

Here  is  the  reflectivity  of  ring  i  and  the  photorefractive  medium  has  length  /. 

We  wish  to  examine  the  stability  of  the  stationary  states  of  this  .system  when  it  is  driven  by  two 
spatio-temporally  distinct  signals  with  similar  input  intensities.  There  are  four  ca.ses  to  consider;  1 ) 
there  is  no  oscillation,  2)  one  of  the  input  signals  oscillates  in  one  ring,  3)  one  input  signal 
oscillates  in  both  rings,  and  4)  one  input  .signal  o.scillates  in  one  ring  and  the  other  input  signal 
oscillates  in  the  other  ring.  The  possibility  of  two  input  signals  oscillating  in  a  single  ring  need  not 
be  considered  since  analysis  of  a  single  ring  resonator  pumped  by  two  input  signals  shows  that  the 
only  stationary  state  corresponds  to  the  stronger  input  signal  oscillating,  with  the  weaker  input 
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signal  completely  suppressed.  (In  the  degenerate  case  where  the  two  inputs  have  the  same  mlen.dt\ 
a  linear  combination  of  both  signals  may  oscillate.) 

The  details  of  the  linear  stability  analysis  are  rather  lengthy  [8].  It  is  found  that  the  only  stable 
solution  is  the  desired  state  w'here  one  signal  oscillates  in  one  ring,  and  the  other  signal  oscillates  in 
the  other  ring.  This  is  a  stable  state,  provided  the  overall  gain  puts  the  circuit  above  threshold,  as 
long  as  the  larger  input  signal  picks  the  ring  with  lower  losses.  Otherwise  the  signals  simply 
switch  rings.  The  degenerate  case  where  the  input  signals  have  equal  intensity  [29).  and  the  two 
rings  have  equal  losses  is  somewhat  arbitrary  since  there  is  a  family  of  .stable  solutions  described 
by  •^i;/-hi  =  =  const.  Numerical  studies  of  the  full  time  dependent  set  of  equations  verify 

that  when  random  noise  is  added  to  the  input  signals  the  resulting  values  of  c  are  uniformly 
distributed. 


III.  Spatio-temporal  dynamics  in  photorefractive  ring  oscillators 


The  self-imaging  ring  resonator  is  a  highly  degenerate  device.  Since  any  transverse  electric  field 
profile,  neglecting  diffraction,  is  imaged  onto  itself  after  one  round  trip  the  eigenmodes  of  the 
empty  cavity  form  a  continuous  set.  !n  this  situation,  where  the  passive  optics  play  only  a  minor 
role  in  determining  the  mode  structure,  the  impact  of  the  nonlinear  elements  comes  to  the  forefront. 
The  self-imaging  ring  is  therefore  an  interesting  test  bed,  both  for  studying  the  effects  of 
nonlinearity,  and  for  leammg  how  to  use  nonlinearity  to  control  the  spatio-temporal  behavior  in  the 
resonator. 

The  single  transverse  mode  ring  resonator  with  photorefractive  gain  and  loss  has  been  shown  to 
exhibit  bistability  and  self-pulsing  [1,2].  The  interesting  dynamical  behavior  in 
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Gain  Plane  Loss  Plane 


Fig.  ni.l  Imaging  ring  resonator  with  photorefractive  gain  and  loss.  All  lenses  are  f=100  mm  with 
a  spacing  of  2f.  The  gain  and  loss  pumps  are  from  a  cw  Argon  laser,  X=5 14  nm. 

this  system  is  due  to  the  competitive  interaction  of  the  gain  and  the  loss.  On  the  other  hand,  in  the 
imaging  resonator  with  a  continuum  of  transverse  modes,  (see  Fig.  III.  1)  the  interaction  of  the  gain 
and  the  loss  assumes  a  cooperative  character  in  the  following  .sense.  Consider  first  what  happens 
when  only  the  gain  pump  is  turned  on.  If  the  gain  pump  has  a  Gaussian  profile  then  the  o.scillating 
mode  will  also  assume  a  Gaussian  profile,  in  the  gain  crystal,  in  order  to  maximize  its  overlap  with 
the  gain  pump.  However,  the  phase  of  the  oscillating  mode  is  unconstrained  since  it  has  no  effect 
on  the  energy  transfer  efficiency  in  the  gain  medium.  The  Fourier  transform  of  a  field  with 
Gaussian  intensity  profile  but  random  phase  is  a  rather  random  looking  pattern.  Now  turn  on  the 
loss  pump  that  interacts  with  this  random  pattern.  When  the  random  pattern  contracts  to  a  localized 
mode  it  maximizes  its  peak  intensity  and  saturates  the  loss  mechanism.  There  is  a  limit  to  how  far 
the  mode  will  contract,  since  if  it  becomes  too  small  its  Fourier  transform  will  become  too  big  and 
no  longer  be  optimally  matched  to  the  gain  pump  in  the  conjugate  plane.  Thus,  if  the  gain  pump  is 
a  Gaussian  with  spot  size  h;,  where  is  the  confocal  mode  size  of  the  passive  cavity, 

and  the  loss  pump  is  somewhat  larger  and  transversely  uniform  then  the  gain  and  the  loss 
interactions  will 
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Fig.  III. 2.  Collapse  of  the  transverse  mode  struetiire.  The  first  frame  shtnvs  the  mode  protiii.' 
immediately  before  the  loss  primp  is  turned  on.  The  following  frames  are  separated  1\\  Feti  ^e*... 
witli  the  gain  plane  on  the  left  and  the  loss  plane  on  the  right. 

cooperate  spatially  tt'  from  an  oscillating  mode  with  Gaussian  spot  si/e  vc  .  bu.  arbitraiy  location 
in  the  plane  of  the  loss  inteniciion.  An  e.xperimental  recording  ot  this  effect  i''  show  n  in  Fig  HI. 2 

r>l. 

However,  once  the  mode  profile  has  collapsed  it  begins  to  wander  about  in  the  transNerve  plane 
1 10].  .An  example  of  this  motion  is  shown  in  Fig.  111..'^.  Depending  on  the  detail''  id  the  CiO  iix 
alignment,  and  the  gain  and  loss  pump  intensities,  the  spot  either  drifts  to  the  edge  of  the  av  ailable 
aperture  and  disappears  before  reappearing  again  in  its  original  position,  or  eKe  it  executes  a 
cyclic  motion  wholly  within  the  aperture.  In  the  latter  case  the  spot  does  not  move  in  a  closed 
circle  but  lather  appears,  moves  a  short  distance,  and  then  'isappears.  before  repeating  its  motion 
It  takes  a  few  seconds  for  the  spot  to  execute  a  single  traverse  across  the  aperture.  The  spot  will 
tvpicallv  choose  a  new  drift  path  after  about  10  cycles  along  any  given  path.  The  wandering 
motion,  alternating  between  paths  that  are  each  repeated  a  few  times,  continues  indefinitely.  The 
spot  motion  appears  to  be  determined  by  the  cavity  alignment  since  small  adjustments  to  one  ot 
the  resonator  mirrors  change  the  spot  trajectory  noticeably. 


[  ly,  111,3.  iApcriniciiial  cibscrvation  ot  spot  The  frames  were  reecnJed  at  0.  3.  anJ 

10  seconds  ot  elapsed  time,  stariine  at  the  upper  left  and  ending  at  the  louer  right. 


The  sptU  motion  has  been  simulated  by  nutiierieal  ealculations  of  the  resiMialor  dsnamies  |21  h 
The  (iptieal  fiel  !  and  the  photorefractive  gratings  are  disereti/ed  on  a  square  lattice.  We  assume  no 
crossing  o!  the  pump  and  signal  beams  in  the  gain  or  k)ss  crystals,  but  retain  the  a\ial  dependence 
of  the  field  and  gratings  in  each  crystal.  Thus  there  is  no  nonlocal  coupling  within  the  crystals 
I'he  transxeise  coupling  which  leads  to  mode  collapse  arises  from  placing  the  crystals  in  spaiialty 
coniugale  planes,  [f  the  ca\  ity  is  assumed  to  be  perfectly  aligned  the  sjm'I  lorms  ui  a  location 
deter  rttmed  bv  the  ma.'vimum  of  the  initial  seedint:.  and  is  stationarc.  Misalienment  id  the  ca\  it\  is 

-  Co  .  . 

modeled  by  intriulucing  a  phase  wedge  in  the  cavity  <I>(  a.  v)  =  e.\p/(c'>,,\  +0  \  ).  w  here  arc 

the  iransxcrsc  coordinates  and  S  .S  characten/e  the  slope  of  the  wedge  Placing  the  wedge  next 
t('  the  gain  crysial  results  in  a  continuous  motion  of  the  spot  parallel  to  the  wedge  giadienl 
Vc|)  ~  b'  \  b  V  .An  analy  tical  study  of  the  spot  molimi  |21 )  lor  the  case  ('f  one  tians\eioc 
dimension  predicts  that  for  sntall  cavity  misalignment  the  drift  rate  is 
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Figure  III. 4  Rate  of  spot  motion  vs.  cavity  misalignment  as  calculated  by  direct  numerical 
simulation.  For  small  S  i  e  rate  depends  linearly  on  6.  The  actual  rate  agrees  with  Eq.flll.l )  to 
within  about  30%,  depending  ^n  the  values  of  the  resonator  parameters. 


f  d 

dt  rln(/^/7,„,)’ 


(Ill.l) 


where  .v  is  the  position  of  the  spot  in  the  loss  plane,  r  is  the  time  constant  of  the  gain  medium.  /,. 
is  the  maximum  pump  beam  intensity  and  is  the  maximum  oscillating  intensity  at  the  entrance 
to  the  gain  medium.  In  deriving  Eq.  (III.l)  the  time  constant  of  the  loss  crystal  has  been  put  to 
zero.  In  the  actual  experiment  the  time  constant  of  the  loss  crystal  w-as  approximately  20  times 
shorter  than  that  of  the  gain  crystal,  so  this  is  a  reasonable  approximation.  Comparison  of  Eq. 
(III.l)  with  the  drift  rate  found  from  direct  numerical  simulations  is  shown  in  Fig.  II1.4.  The  drift 
rate  is  indeed  linearly  proportional  to  d  for  small  5.  it  is  surprising  how  small  8  must  be  before  a 
linear  dependence  of  drift  rate  on  cavity  misalignment  is  obtained.  This  indicates  that  a  higher 
order  theory,  which  allows  for  a  nonlinear  dependence  on  5  may  be  necessary  to  get  good 
agreement  between  the  analytical  and  numerical  descriptions. 

It  appears  clear  from  the  analytical  and  numerical  studies  that  the  motion  of  the  localized  mode  is  a 
direct  consequence  of  cavity  misalignment.  Even  if  the  misalignment  is  made  very  small  the  spot 
will  continue  to  move,  albeit  slowly.  In  order  to  u.se  this  type  of  resonator  for  implementation  of  a 
neuromorphic  optical  processor  we  wish  to  eliminate  the  spot  motion.  We  have  very  recently  found 
a  configuration  involving  a  resonator  with  three  photorefractive  crystals  that  eliminates  the  snot 
motion  in  our  computer  simulations.  On  going  work  will  concentrate  on  experimental  vtriiication 
and  analytical  study  of  this  new  configuration. 
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IV.  Transverse  instability  of  counterpropagating  beams 


Transverse  modulational  insiabilily  cd' two  counterpropagating  beams  m  nonlinear  optical  media 
has  been  observed  in  atomic  vapors  {13-151,  liquid  crystals  [16.17],  and  quite  reeentl>  m 
photorefracti\ e  crystals  1 18],  The  instability  manilests  itself  in  the  far  field  b\  the  appearance  ol 
conical  rings,  pairs  of  sp,its  or  a  he.xagonal  array  of  spots  about  the  counterpropagating  beams. 
The  corresponding  intensity  modulation  in  the  near  field  has  the  characteristic  spatial  scale 
/  =  2;r/0  A,,,  where  0  is  the  small  angle  between  the  generated  satellite  beams  and  the  primar> 
beams.  The  angle  0  at  which  the  threshold  of  the  absolute  instability  is  a  minimum  may  be  tiuind 
for  all  of  the  above  materials  by  a  linearized  analysis  of  the  equations  of  motion  [  12.1^.20.22-26] 
On  the  other  hand,  the  nature  of  the  nonlinear  stage  of  the  instability,  the  details  of  the  resulting 
patterns,  and  their  spatial  and  temporal  stability  appear  to  vary  u  idely  among  different  materials  In 
atomic  vapors,  cones,  pairs  of  spots  and  hexagons  may  all  be  observed  depending  on  the 
frequency  and  intensity  of  the  primary  beams.  In  liquid  crystals  the  instability  leads  to  a  stable 
hexagonal  pattern,  w  hereas  the  recent  observations  in  photorefractive  media  indicate  a  somewh.ii 
different  behavior.  Here  the  hexagonal  pattern  is  well  defined,  hut  temporally  and  spatially 
unstable,  such  that  it  appears  as  a  cone  in  a  time  averaged  recording.  The  cones  observed  in  atomic 
vapors  may  also  be  due  to  unstable  hexagons  although  no  time  resolved  data  is  presently  available. 
There  is  some  reason  to  expect  this  since,  as  was  shown  in  [  15].  the  cone  may  he  collapsed  into  a 
hexagon  by  injection  of  a  weak  beam  at  0  , 

Although  the  general  nature  of  the  instability  in  photorefractive  media  is  similar  to  the  well  known 
Kerr  nonlinearity  case,  the  details  of  the  instability  threshold  conditions  differ.  This  is  due  ti' 
several  differences  between  the  photoref:  ctive  and  Kerr-type  nonlmearities.  In  Kerr-type  media 
the  nonlinear  part  of  the  refractive  index  n  is  not  (or  is  only  weakly  )  dependent  on  the  angle  0 
between  the  interacting  waves  (sec  Fig.  I\  .  1 1  or.  equivalently  ,  on  the  wavevector  k  =  0k.  of  the 
grating  written  by  the  pumping  waves  w  ith  their  sidebands  ( w  ith  and  B..  w  ith  /T.  i  and  so 
the  characteristic  angle  0  is  the  result  of  interplay  between  the  Kerr  nonlinearity  and  diffraction. 
The  nonlinear  coupling  coefficient  (analog  of  tu  i  in  phomrefractive  media  is  strongly  dependent  on 
the  value  of  k  and  so  material  properties  of  the  crystal  must  come  into  play,  imposing  their  viwn 
characteristic  spatial  scales.  Furthermore,  the  Kerr  nonlinearity  corresponds  ti^  a 


Nonlincur 

Medium 


(n) 


(b)  (c) 

Figure  IV.  1,  Transverse  instability  of  counterpropagating  beams,  a)  Geometry  of  the  optical 
interaction;  incident  pump  beams  while  (F'.pfi.,)  are  spontaneously  generated 

satellites  b)  transmission  gratings  couple  the  beam  pairs  )}.  and 

)}  ,  and  c)  reflection  gratings  couple  the  beam  pairs  )}. 

nonlinear  change  of  pha.se  ( ;i,  is  purely  real)  whereas  the  photorefractive  nonlinearity  is  in  general 
complex  corre.sponding  to  both  amplitude  and  pha.se  changes.  The  magnitude  of  the  real  part  of  the 
coupling  coefficient  can  be  enhanced  by  applying  an  external  electric  field  to  the  photorcfracti%c 
crystal,  but  in  general  it  is  impossible  to  completely  eliminate  the  imaginary  part.  Also  ihc 
dependence  of  the  real  and  imaginary  parts  on  the  value  of  k  is  different.  In  addition 
photorefractive  media  are  characterized  by  strong  amplified  incoherent  scattering  (fanning  i  which 
leads  to  a  dependence  of  the  amplitudes  of  the  primary  beams  on  the  axial  coordinate,  even  in  the 
linear  stage  of  the  instability. 

The  details  of  the  analysis  in  photorefractive  media  depend  on  whether  the  in.stability  i>  due  to  the 
formation  of  transmission  or  reflection  gratings.  The  interacting  waves  arc  written  in  the  form 

Fir  ,zj)=  F,j(z)[]  +  F^^txpiik  r  -  iQl)+  F  ,cxp{-ik  r  +  iQ  !)]cxp{ik^^z- 

I  "  (IV.  i ) 

B(T\,z.t)  =  B„(c)Il  +  iB,|  exp(/^^r  -  iXJt )+  B  ^exp{-ik  r  +  iQ‘ ^)]e\p^-ik^~  -  /m,/) 
which  leads  to  the  equations  of  motion 
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{dj^^  +  /AJF  ,  =  ,  +  fi’,  )]• 

(^/*-/A-,)f;,  =-/-i^[f ,  +  f;,  +^(fl ,  +  f‘,)]. 
\  + 

[djdz  -  ik,)B,,  =  +  /"-',  +  <7(fi  ,  +  B’, )]. 

[d/Bz  +  ik,)B[,  =  /-^[F  ,  +  r,  +  ^(fi  ,  +  f:,  )] 

1  +  <y 


for  transmission  gratings  and 

iLl^iyS^F 


dB^ 

d:. 


=  -iy 


\  +  q 


{d/az + ik,)F^,  = '  ]• 

{<?/*- /A;, )f:,  =-/y^[(7-'  -7’)f:,  +y’F:,J. 

{dldz-ik,)B,,  =-iz^[yF^,+Y-'F:,+{y-’-y)B,\ 
{dj^z  +  /A-Jf;,  =  '■  j^[7'  ^.1  +  yB'.^  +(/'  -  7)5',]. 


(IV.3) 


for  reflection  gratings.  In  Eqs.  (IV. 2)  ?(;:)  =  |Fo(c)//\)(c)|'.  Aj  =  A; /2A(,.  and  y{k  .Q)  is  the 
material  and  frequency  shift  dependent  complex  coupling  coefficient.  In  (IV. 3)  y  =  7(2A„.Q  =  0) 
and  r=yn+^T^ . 


The  dispersion  relationship  governing  the  instability  threshold  is  found  by  solving  Eqs.  (1V.2. 
IV. 3)  with  the  boundary  conditions  F^,(c  =  0)  =  F_,(0)  =  0.  B,,(z  =  !)=  B  ,(1)  =  0.  In  \hc  case  of 
transmission  gratings  Eqs. (IV. 2)  may  be  solved  in  closed  form  for  arbitrary  axial  variation  of 
q{z)  We  find  that  the  dispersion  relation  is  given  by 
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(1-A)(1-^:)--4A=0, 

I  l  +  qiz)^in{kj) 

A  =  (- Aj/OJ dz  -  sinh(.v:). 

I  l  +  q(z.)  sin{;,/) 

A,,  =  (2}<-,/.v)J^/ 


f,,,  ■?(-^)  ^i"l^f-/>l,i„h[,(,-,)]. 

l+^(:)  sm{kj) 


\  =  {2-ikJs)\dz 


1  sin(A:,,c)  , 


l  +  ^(.-}  sin(y) 


sinh(.vc). 


(iV,4l 


For  ^  const,  gqj.  (|v.4)  reduces  to 

(<7  +  r/'')  +  |(.v/A',)  '  -  {sjkj)  3in(A:j/)sinh(,v/)  + 2cos(A:/)cosh(.v/)  =  0. 


(1V.5) 


If  we  impose  the  restriction  of  purely  real  y  and  Qi  =  0  this  equation  corresponds  to  the  know  n 
results  for  Kerr  media  [12.22-26], 


Several  of  the  lowest  branches  found  from  numerical  solution  of  Eq.  (IV.5)  are  shown  in  Fig. 
IV. 2.  Taking  q  =  const,  is  usually  not  a  good  approximation  in  photorefractive  media  due  to  the 
influence  of  passive  absorption,  and  beam  fanning.  The  amount  of  beam  fanning  may  be  reduced 
by  w'orking  with  narrow'  beams,  although  in  the  transmission  grating  geometry  the  beam  diameter 
cannot  be  made  too  small  due  to  the  requirement  of  having  at  least  a  few  modulation  fringes  across 
the  main  beams.  Numerical  solution  of  Eqs.  (IV. 4)  with  ^(r)  =  where  cr,„  =  9cm  ' 

accounts  for  passive  absorption  and  fanning  as  measured  with  mm  sized  beams  in  a  piece  of 
BaTiOa,  shows  that  the 
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Fig.  IV. 2.  Dispersion  curves  for  q(z)=const.:  a)  focusing  branch  with  Re(*/)>0  and  b) 
defocusing  branch  with  Re(7)  <  0.  The  curves  are  labeled  as  i)  q=l.  ii)  q=10.  and  iii)  q=100. 

-  M- - arg(7), 

necessary  coupling  constant  for  observation  of  instabilities  is  increased  by  about  a  factor  of 


The  recent  experimental  observation  of  the  modulational  instability  in  photorefractive  media[18] 
used  a  reflection  grating  geometry.  This  is  advantageous  since  it  is  then  possible  to  use  very 
narrow  beams,  which  reduces  fanning  tosses.  Eqs.  (IV. 3)  describing  the  reflection  grating  case  are 
more  complex  than  those  describing  transmission  gratings  and  general  solutions  for  q  =  q{z)  are 
not  available.  However,  by  choosing  the  boundary  conditions  such  that  =  |fio(0i  I'esults  in 

q-\  throughout  the  photorefractive  medium  since  Eqs.  (IV. 3)  have  the  first  integral 
=  const,  for  any  value  of  y.  Putting  q  =  \  in  Eqs.  (IV. 3)  results  in  the 
dispersion  relationship 
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